Internal time, test clocks and singularity resolution in dust-filled quantum cosmology 
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The problem of time evolution in quantum cosmology is studied in the context of a dust-filled, 
spatially flat Friedmann-Robertson- Walker universe. In this model, two versions of the commonly- 
adopted notion of internal time can be implemented in the same quantization, and are found to yield 
contradictory views of the same quantum state: with one choice, the big-bang singularity appears 
to be resolved, but with another choice it does not. This and other considerations lead to the 
conclusion that the notion of internal time as it is usually implemented has no satisfactory physical 
interpretation. A recently proposed variant of the relational-time construction, using a test clock 
that is regarded as internal to a specific observer, appears to provide an improved account of time 
evolution relative to the proper time that elapses along the observer's worldline. This construction 
permits the derivation of consistent joint probability densities for observable quantities, which can 
be viewed either as evolving with proper time or as describing correlations in a timeless manner. 
Section |T] reveals whether the singularity is resolved or not. 

PACS numbers: 98.80.Qc, 04.20.Cv, 04.60.Ds 



I. INTRODUCTION 

It has long been appreciated that time evolution in 
generally-covariant theories such as general relativity is, 
on the face of it, a gauge transformation and, contrary to 
everyday experience, should therefore be unobservable. 
Reviews of this 'problem of time' are given, for example, 
in [ll-Q , and textbook discussions may be found in [1, Q . 
In recent years, detailed studies of quantized models of 
cosmology (see, e.g. [rj|-[i"o| for reviews) have demanded a 
practical solution to this problem, and a certain notion of 
'internal time' introduced by Rovelli [lll - fl3j and further 
developed in, for example, [Mrfl7| has been quite widely 
adopted. The 'evolving constant of the motion' construc- 
tion advocated in these papers is a particular implemen- 
tation of the idea of relational time, according to which 
time evolution in covariant theories can be described only 
relative to the values assumed by some physical quantity 
that is chosen to serve as a clock. In simple cosmolog- 
ical models, the clock is typically a scalar field, say </>, 
and in the quantum theory one obtains a wavefunction 
ij)(v;ip), where v is a variable describing the geometry of 
the model universe, as the solution of a Schrodinger-like 
equation in which the parameter (p plays the role of time. 
In this paper, we take v to be the volume of a spatially 
compact universe, in which case \if)(v; <p)\ 2 is interpreted 
as a probability density for the volume at the time when 
the scalar field <fi takes the value cp. 

Recently [l8|, we have argued that this interpretation 
is hard to sustain because, for the reasons summarized 
in section IIV Al below, the values assigned to (p cannot 
be regarded as readings obtained by inspection of the 
'clock' (f>. Indeed, according to the usual rules of quan- 
tum mechanics, the scalar field is not an observable at 
all, which is especially disconcerting if it constitutes the 



entire matter content of the universe. In view of this dif- 
ficulty, we proposed a variant of the relational-time idea, 
which makes use of a small test clock, regarded as internal 
to some observer from whose point of view the universe 
is to be described. In this alternative construction, the 
test clock is unobservable, and we obtain a wavefunc- 
tiorQ 4>(v,(f>]T), evolving with a time parameter t which 
is not a clock reading, but corresponds classically to the 
geometrical proper time that elapses along the observer's 
worldline. This new wavefunction yields a joint proba- 
bility density, evolving in textbook fashion with proper 
time r, for v and which are now genuine observable 
quantities. 

In this paper, we first wish to investigate whether, 
despite the above-mentioned difficulties, the internal- 
time wavefunction ip(v;ip) can be interpreted as express- 
ing a correlation between two quantities (the volume 
and scalar field), both of which are observable in some 
suitably broadened sense. Specifically, we ask whether 
| ip( v i f) 1 2 can be regarded as a conditional probability 
density for the volume, given that the scalar field has 
been determined to have the value (p. To that end, we 
study a simple model of an homogeneous universe filled 
with pressureless matter, described, following Brown and 
Kuchaf [l9j, by a single scalar field. The classical version 
of this model is introduced in section [TTl and we find 
that two complementary notions of internal time can be 
straightforwardly defined, using either the volume or the 
scalar field as a clock. These two internal times carry over 
to the quantized theory, as discussed in section ITO1 and 
we use them in section HVl to compute two corresponding 
probability densities. // the conditional-probability in- 
terpretation is feasible, then these ought to represent, on 
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1 Later on, we will adjust the notation in which various wavefunc- 
tions are expressed, so as to maintain some important distinc- 
tions, which will be made precise in due course. 
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the one hand the conditional probability density for the 
volume given some value of the scalar field and, on the 
other hand, the conditional probability for the scalar field 
given some value of the volume. We find, however, that 
these probabilities are inconsistent. In fact, the two no- 
tions of internal time lead to two mutually contradictory 
views of time evolution: if the volume is used as inter- 
nal time, the quantum theory appears to reproduce the 
classical big-bang singularity, whereas, if the scalar field 
is used as internal time, exactly the same quantum state 
appears to describe a universe in which the singularity is 
replaced by a bounce. This result serves to strengthen 
our conviction that the internal-time formalism, while 
resulting from perfectly sound mathematical manipula- 
tions, has no satisfactory physical interpretation. 

In section |V1 we augment the model by the addi- 
tion of a test clock. As already indicated, the modi- 
fied relational-time construction proposed in (Tg] leads 
to a joint probability density for the volume and scalar 
field, both of which which are now bona fide observables, 
which evolves, according to a standard Schrodinger equa- 
tion, with proper time r. Like the usual time parameter 
in non-relativistic physics, r is not itself observable. In 
principle, no set of experimental results records the evolu- 
tion of an observed quantity relative to the time param- 
eter that appears in Newton's laws or in Schrodinger's 
equation. Rather, what is recorded is an observed corre- 
lation between the object of interest and the readings of 
some time-keeping device, and the rationale underlying 
the internal-time formalism is that taking this practical 
fact seriously is the key to solving the problem of time. 
Consequently, it is of some interest to recover a 'timeless' 
interpretation, by supposing that the observer to whom 
our unobserved test clock is internal also possesses a time- 
keeping instrument that is observed. This issue is taken 
up in section IV CI For the model at hand, though not in 
general, it turns out that the role of a time-keeping in- 
strument can conveniently be assigned to the scalar field 
(but not to the volume). When that is done, we obtain 
a genuine conditional probability density for the volume, 
which coincides, to a good approximation, with the func- 
tion \ip(v; tp)\ 2 found from the internal-time formalism 
with ip as the internal time. In that indirect and approx- 
imate sense, the Brown-Kuchaf scalar field emerges as a 
preferred internal time, indicating a bouncing universe. 



II. CLASSICAL DUST-FILLED COSMOLOGY 

A. The model 

As described in detail by Brown and Kuchaf (and 
in earlier work by Brown [20]), an effective hydrodynamic 
description of pressureless matter, or dust, is furnished by 
a collection of scalar fields {</>, Wk, Z k , M} (k = 1,2,3), 
with a Lagrangian density of the form 

£d = -\sf~gM [gTUpUv + 1) , (2.1) 



where the 4-velocity field is = -d^fy + Wkd^Z k . 
When the equations of motion are satisfied, the inte- 
gral curves of U can be interpreted as the worldlines 
of dust particles, and the scalar field cf> is linear in the 
geometrical proper time t that elapses along these world- 
lines. We consider the case of an homogeneous, spatially 
flat Fricdmann-Robertson- Walker universe, with metric 
g = diag(— N 2 , a 2 , a 2 , a 2 ). To be concrete, we take spa- 
tial sections of this universe to be compact, with coordi- 
nate volume J d 3 x = 1. Homogeneity implies that spatial 
derivatives of (j> and Z k vanish, and variation with respect 
to Wk leads also to doZ k = 0. In this special case, there- 
fore, the dust is modeled by a single scalar field (f>, along 
with the Lagrange multiplier M: 

C u = §Na 3 M [N- 2 (d s <f>) 2 - 1] , (2.2) 

where s is an arbitrary time coordinate. The momen- 
tum conjugate to <j> is = N~ 1 a 3 Md s <j), and varia- 
tion with respect to M yields the second-class constraint 
p 2 — (a?M) 2 = 0. Up to a sign, this constraint is trivially 
solved for M, and we can construct the unconstrained 
dust Hamiltonian 

H D = Np^. (2.3) 

Evidently, p^ is the total energy content of the dust, and 
we resolve the sign ambiguity by requiring this to be non- 
negative. 

In the standard way, the Einstein- Hilbcrt action leads 
to a gravitational Hamiltonian, which we write in terms 
of the volume v = a 3 and its conjugate momentum p v — 
-(UnG^N^dsv/v as 

tfgrav = -N{6nG)vp 2 v =: 7VC grav , (2.4) 

G being the usual gravitational constant. The function 
Cg ra v defined by this equation is the gravitational contri- 
bution to the Hamiltonian constraint. 

Taking the lapse function N(s) to be a strictly posi- 
tive, but otherwise arbitrary function, we introduce the 
invariant proper time 

t(s) = f S N(s')ds' (2.5) 
Jo 

and express the Hamilton equations of motion as 

i> = N- l dH /d P4> = 1 (2.6) 

v = N~ 1 dH /dp v = -(12itG)v Pv (2.7) 

Pv = -N^dHo/dv = (6t:G)p 2 v , (2.8) 

where the overdot denotes differentiation with respect 
to t. The total Hamiltonian here is Ho = if grav + Hd 
(distinguished by its subscript from the Hamiltonian of 
an extended model to be considered later) and the energy 
Pcf, is a constant of the motion. Finally, the Hamiltonian 
constraint 

Co :=dH /dN = -(6irG)vpl+ Ptl> = (2.9) 
reproduces the Friedmann equation. 
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B. Dirac observables and internal time 

The equations of motion (|2.6p - (|2.8[) are easily solved. 
Denoting by (f>(v,p v , <f>,p^; i), etc. the phase-space trajec- 
tory that passes through (v,p v , 4>,p<p) at t = 0, we have 

^(v,p v ,(f),p^;t) = <f> + t (2-10) 
v(v,p v ,<j),p$;t) = v(l - 6irGp v t) 2 (2.11) 
p v {v,p v ^(j),p^;t) = p v (l - QTrGpyt)' 1 . (2.12) 

Heuristically, given initial conditions that satisfy the con- 
straint (|2.9[) . these solutions provide a complete descrip- 
tion of the evolution of this simple universe relative to 
the proper time t that elapses along the worldline of a 
comoving observer. Classically, this is possible because 
the equations of motion were derived before imposing 
the constraint. In a quantum-mechanical treatment, this 
time evolution cannot be reproduced, because the Hamil- 
tonian Cq that is supposed to generate it vanishes when 
acting on a physically allowed state. The same difficulty 
arises classically in a more formal treatment of the con- 
strained Hamiltonian dynamics. Here, one has to recog- 
nize that, while the proper time defined in (I2.5[) is invari- 
ant under reparametrizations of the coordinate s with 
N'(s') — N(s)ds/ds' (the remnant, in this symmetry- 
reduced model, of the general coordinate invariance of 
general relativity), it is not invariant under an arbitrary 
change in the undetermined lapse function N(s) unless 
this is accompanied by a compensating reparametriza- 
tion. From this perspective, a change in N(s), and hence 
evolution with respect to t, is a gauge transformation 
generated by the constraint function Co- One has the 
option of fixing a gauge, by specifying once and for all a 
definite function N(s). Classically, at least, the content 
of (|2.10[) - (|2.12p is independent of the actual function cho- 
sen, and can be regarded as physically meaningful. But 
in a manifestly gauge-independent approach, which we 
follow here, genuine physical information is carried only 
by gauge-invariant (Dirac) observables, which commute, 
in the sense of Poisson brackets, with Cq, and are there- 
fore constants of the motion. 

A construction due to RovellifTll [l~3| allows us to ob- 
tain a 1-parameter family of gauge-invariant quantities — 
an 'evolving constant of the motion' — as follows. Denote 
by t v the time at which cf> in (|2.10p has the value ip. Then 
the quantity 

V(ip) := v(t v ) = v[l - ^Gp v {p - 0)] 2 , (2.13) 

in which we suppress the dependence on the phase-space 
coordinates (v,p Vl 4>,p<p), can be interpreted classically as 
'the volume at the time when the scalar field has the value 
(p\ It is easy to check explicitly that {V((p), Co} = for 
each value of the parameter ip. 

A parameter such as ip is commonly referred to in the 
literature as a 'relational', 'emergent' or 'internal' time, 
or simply as 'time', the idea being that the scalar field 
serves as a physical clock, and the function V(ip) de- 
scribes the evolution of the volume with respect to the 



readings of this clock. In the present case, this may seem 
especially apt, in view of the linear dependence (|2.10l) of 
the scalar field on t. In [l9| . indeed, this idea is enshrined 
in the notation: these authors use T for the scalar field 
that we denote by (f>. For reasons that will become clear 
later, however, we wish to maintain a clear distinction be- 
tween the geometrical proper time defined by (|2.5[) and a 
scalar field <p whose equation of motion happens to have 
the solution (12 . 10[) . 

We now wish to define a second family of Dirac observ- 
ables, taking the volume, rather than the scalar field, 
as an internal time. This entails solving (|2.11j) for the 
time t v at which the volume has the value v, and re- 
quires a little care with regard to signs. First, we take 
the volume always to be positive, in contrast to the loop- 
quantum-gravity-inspired treatment described in [2ll l22j | , 
where the physical volume is the absolute value of a 
more fundamental variable, whose sign reflects the ori- 
entation of a co-triad. Next, it follows from the equation 
of motion (|2.8p that the evolution preserves the sign of 
p v , so the quantity a that we provisionally identify as 
a = — sgn(p„)> is a constant of the motion. Inspection 
of the solutions (|2.1ip and (I2.12p then reveals that tra- 
jectories on which p v is negative correspond to an ex- 
panding universe, starting from an initial singularity at 
t = — (GirGlpvl) -1 , while those on which p v is positive 
correspond to a contracting universe and terminate at a 
final singularity at t — (6irGp v ) . Consequently, we can 
take the square root of ([2. lip , with v = v, to obtain 

u 1 ' 2 = v 1 ' 2 - ^Gv x ' 2 p v t Vl (2.14) 

the sign of the square root being determined unambigu- 
ously by the requirement that v is an increasing function 
of t v when p v is negative. Thus, the Dirac observable 
that represents 'the value of the scalar field when the 
volume is x/ is 

:= 0(i„) = <\> + (GvrGp,,)- 1 - {6nGv 1 / 2 Pv y 1 v 1 ' 2 . 

(2.15) 

Again, one may check that {$(^), Co} = for every value 
of v. 

Finally, with a view to quantization, we define the 
Dirac observables 

V := V(0) = v + l2nGvp v 4> + {<6irG) 2 vp 2 v (j) 2 (2.16) 

Y := — (127rG) _1 y'(0) = vp v + 6nGvp 2 v (/) (2.17) 
$:=$(O) = + (6vrGp„)- 1 , (2.18) 

and we note that C grav , defined in (I2.4p . is a constant of 
the motion, and hence also a Dirac observable. We have 
the Poisson-bracket relations 

{V, Y} — V, {V, C grav } = - 12nGY, 
{Y, G grav } = G grav , {$, C grav } - -1 (2.19) 

and the evolving observables can be expressed as 

V{<p) = V — UirGYip - 67rGC grav ^ 2 (2.20) 
= $ + a{-6TTGC grav y 1/2 v 1/2 , (2.21) 
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where we now identify 



-sg 



n(v^ 2 p v ). 



(2.22 



The factor v 1 / 2 here, which follows from (|2.15p appears 
inessential, but it will prove convenient to retain it. 

We now have two notions of evolution with respect to 
internal time, and it will be useful to identify the gener- 
ators of these evolutions. Indeed, we easily discover from 
(l2~T9l) - (l2T2T1) that 



dip 
d$(z/) 
di/V2 



{V(<p),H v } 



with 



Hip . — Cgrav 



H v := ai-C^/^G) 1 / 2 . 



(2.23) 
(2.24) 

(2.25) 
(2.26) 



On any classical trajectory, <&(v) is just the inverse of 
V(ip), but we shall see that this inverse relationship is 
not preserved in the quantum theory. 



III. QUANTUM DUST-FILLED COSMOLOGY 

We consider a quantization scheme of the Wheeler-de 
Witt type in which, in the first instance, the canonical co- 
ordinates (v,p v ,<f),P(f,) are promoted to operators acting 
in an auxiliary (or kinematical) vector space. A conve- 
nient representation is that in which the operators v and 
Prj, act by multiplication on wavefunctions ^(v, en), while 
their conjugate variables act by differentiation: 



p^iv^eu) = e D *(w,e D ) 
0*(u,e D ) = iW en ^(v,e ). 



(3.1) 



The notation cd reflects the fact that corresponds to 



the energy content of the dust. As in [18j we follow the 
authors of [U [22[ in choosing for the gravitational con- 
straint the operator ordering 



-6irGp v vp v . 



Then the Hamiltonian constraint equation 

[(4/\ 2 )d v vd v +e D ]y(v,e D ) = 0, 

where 



A := 



3irGh 2 



1/2 



With the definition 



z := 



Aer/V/ 2 , 



(3.2) 
reads 
(3.3) 

(3.4) 

(3.5) 



the general solution to this equation can be expressed as 



(3.6) 



where, in terms of the usual Hankcl functions, we write 
H + (z) = -iH (2 \z) and H_(z) = iH^\z). 

The classical phase space consists of two disjoint re- 
gions, distinguished by the discrete variable (|2.22p . con- 
taining expanding (a = 1) and contracting (a = — 1) 
trajectories. (We ignore, for now, the hyperplane p v = 0, 
on which the volume is constant, but see the comment 
following (14.14(h ) The corresponding operator is 



sgn 



ihXe^ 2 d 
2 dz 



and we see from the integral representation 



2 f c 
H ± {z) = - 

t Jo 



d£ exp (=Fi,zcosh£) 



(3.7) 



(3.8) 



that H + (z) is a linear superposition of eigenfunctions of 
id z with positive eigenvalues, while H-(z) is a superpo- 
sition of ei genfunc tions with negative eigenvalues. More- 
over, since H + (z) = H-(z), these two components of ^ 
are orthogonal if, for example, we choose the inner prod- 
uct 



^ ^ 

de D / dv^i(v,e B )d v ^ 2 (v,ei))- 
h 



(3.9) 

While we will not make direct use of this inner product, 
these consideration provide the heuristic motivation for 
our actual choice of the physical Hilbert space "Hphys- 
Clearly, a solution of the constraint equation is specified 
by a pair of functions V^^d), and we take H p h ys to be 
the direct sum of two copies of L 2 (R + ,deo), with the 
inner product 

/*°° 

(V'l^Ws = / de D \ip + (e u )\ 2 + \iIJ-{£b)\ 2 
Jo 1 

(3-10) 

Then the operator a, which is a constant of the motion, 
has a well-defined action in "Hphys, namely 



-1>- 



(3.11) 



With the exception of H u , defined in (|2.26p , the observ- 
ables we need to consider act independently in the two 
subspaces, and we will use the shorthand Oip(ei)) to rep- 
resent an action of the form 



U-( £ d)/ W-(e D ) 



(3.12) 



The operators v, p v and (j> defined in (|3.1I) do not have 
any well-defined action in "Hphys; that is, they do not 
act on a linear combination of H+{z) and H_(z) with 
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z-independent coefficients to produce another such func- 
tion. As one might expect, however, we can convert the 
classical Dirac observables V, Y and (with slightly more 
difficulty) $ given in (|2.16l) - (|2.18|) into operators that do 
act in 'Hphys- For V and Y we choose the operator order- 
ings 



V = v + 6irG(vp v + p v ii)4> 

Y := \{vp v +p v v) - C grav (, 



av V 



(3.13) 
(3.14) 



with Cgrav given by (|3.2|) . These three operators have 
the commutation relations [A, B] = ih{A, B}, with the 
Poisson brackets shown in (|2.19|) , and it follows that the 
Heisenberg equation of motion (12.23)) is promoted to 



dip 



(3.15) 



with H v = (7 grav and V(ip) the operator version of (|2.20l) . 
With the use of Bessel's equation (zd z zd z + z 2 )H±(z) = 
0, it is straightforward to find the action of these opera- 
tors in 'Hphys : 



4 d d 

V den oen 
d 



YiP(e B )=ihe 1 I / 2 ^-e 1 ^(ev) 
H v ip(e-D) = C SIa , v ip(ex>) = -er>i>(er>)- 



(3.16) 

(3.17) 
(3.18) 



With the inner product (|3.10|) , V and Y are clearly sym- 
metric, and H v , which acts by multiplication, is self- 
adjoint. 

Construction of an operator $ corresponding to (|2,18[) 
is not quite straightforward, because p v = —ihd v does not 
have a well-defined inverse. Consider, however, a wave- 
function 9(v, £d) = i/j(€d)Co{z), where Co is any Bessel 

function of order 0. We define an operator p^ 1 by writing 
_d${v,e B ) ih\ 2 



$*(v,e D ) = ih- 



in— Co(z) 



ih 
2^> 



V(v',en)dv' 



V , (cd)A(z), 



where 



A{z) = zd z C (z) + I C (z')z'dz' 



zq is a constant, and vo — z 2 /X 2 e^. By virtue of Bessel's 
equation, we have 



A(z) - zd z C (z) - / d z , [z'dMz')] dz' 

J z 

— z C (zo), 

and this vanishes, provided that we choose zq to be a 
zero of zC' (z), which is at a complex infinity in the case 
of the Hankel functions. In this way, we obtain 



This operator is symmetric under the inner product 
(|3.10[) . and has the commutator [&,C grav } = —ih, 
in agreement with the last Poisson-bracket relation of 
(|2.19[) . Consequently, the equation of motion (|2.24[) be- 
comes the operator equation 



d^/2 



(3.20) 



where 



/ 2e \ 1/2 

H v = aG Vl G^(e D )=f^|J ^(cd) (3.21) 



$-0(e D ) = ihd eD ip(e u ). 



(3.19) 



is self-adjoint, since it acts independently by multiplica- 
tion in each subspace. 



IV. INTERPRETATION OF INTERNAL TIME 
A. Conditional and joint probabilities 

We argued in 18] that, while an operator such as V(ip) 
is a perfectly good Dirac observable, the internal time pa- 
rameter ip cannot bear the interpretation that one would 
like to place on it (and which indeed is placed on it, in a 
significant part of the literature). 

At the level of the classical equations of motion (12. 6[) - 
(I2.8[) and their solutions (I2.10l) - (|2.12l) . it seems quite fea- 
sible to say that V(ip) is the volume 'when' the scalar 
field has the value ip, provided that one does not enquire 
too closely about the instant of time 'when' this pair of 
values is realized. But even classically, if one follows the 
systematic procedure of constructing a reduced, physi- 
cal phase space, whose points are gauge orbits in the 
constraint manifold, one finds that this physical phase 
space is 2-dimensional. The corresponding configuration 
space is 1-dimensional: there do not exist two indepen- 
dent physical quantities which might simultaneously be 
determined to have the values ip and V(ip). 

Quantum-mechanically, the configuration space on 
which states in "Hphys are defined is again 1-dimensional. 
The parameter ip, which labels the family of Dirac ob- 
servables V(ip), cannot be construed as a value obtained 
by observation of a physical clock (the scalar field), be- 
cause there is no operator acting in H p h ys , independent 
of V, to represent any such observable clock. In par- 
ticular, the operator $ cannot serve serve this purpose, 
for two related reasons. First, $ does not commute with 
V(ip) for any value of ip, so the rules of quantum me- 
chanics do not allow simultaneous measurements of the 
quantities represented by these two operators. Second, 
the reason these operators do not commute is that they 
were constructed through the 'evolving constant of the 
motion' algorithm, and are thus quite different from the 
kinematical operators v and </>. Classically, to interpret 
ip as a value of $ is to interpret V(ip) as "the volume at 
the time when the value ip is assumed by 'the scalar field 
at the time when the volume is zero' " . This is, of course, 
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incoherent, and is made no less so by the transition to 
quantum mechanics. Similar remarks apply, of course 
to the parameter v that labels the family of observables 
$(y): it cannot be construed as the result of a measure- 
ment of the volume. Given that these parameters cannot 
be construed as values obtained from measurements, it is 
hard to see that they have any physical meaning at all. 

We now wish to strengthen this conclusion by con- 
sidering the possibility that, notwithstanding the argu- 
ments just given, the internal-time formalism might be 
construed as yielding a joint probability distribution that 
represents a correlation between observable quantities, in 
this case the volume and the scalar field. That is to say, 
we will try to extend the notion of 'observables' by drop- 
ping the requirement that they be represented by mutu- 
ally commuting operators in "H p hys- This is, in particular, 
a timeless interpretation, in which one might decide to 
ignore conundrums concerning the times at which spe- 
cific values of the observables are realized. We will show, 
though, that in general no such probability distribution 
exists. 

The idea is this. If V(tp) is regarded as the Heisenberg- 
picture operator associated with evolution in an internal 
time tp, generated by the Hamiltonian H v , then we can 
construct the corresponding Schrodinger-picture wave- 
function 



i\)(y, tp) = exp{-\H ip tp/h)tp{v,Q), 



(4.1) 



where i/j(v, 0) is a suitable transform of ^(cd) on which 
V (that is, V"(0)) acts by multiplication: Vip(D,0) — 
v>vl>(v,0). According to the usual interpretation, the ob- 
ject 



V{v;tp) := \iP + (v,tp)\ 2 + \ip_(v,tp)f 



(4.2) 



is the time-dependent probability density for obtaining 
the value v from a measurement of the volume performed 
at 'time' tp. In particular, this probability has the time- 
independent normalization V(v; tp)dD = 1. Accord- 
ing to the foregoing discussion, the problem with this is 
that we have no idea what is meant by 'performing the 
measurement at time tp\ 

The problem might be alleviated if we could reinterpret 
(|4.2p as a conditional probability density V(v\tp) for the 
volume, given that some other quantity (which we hope 
to identify with the scalar field) has been determined to 
have the value tp, even though that other quantity does 
not appear explicitly in our formalism. An interpretation 
of that kind requires the existence of a joint probability 
density V(y, tp), with the normalization 



dv / dtpV(u, tp) = 1, 



such that 



V{D\tp) 



V(u,tp) 



^V{v',tp)dv' 



(4.3) 



(4.4) 



If this joint probability density does exist, we can use it 
to find the conditional probability density 

V(v,<p) 



2{tp\v) = 



(4.5) 



for the scalar field, given that the volume has been de- 
termined to have the value v. 

For the model under consideration, it happens that we 
already have such a probability density to hand, namely 

V{tp\v) :=\f + {tp,v)\ 2 + \f_{tp,v)\ 2 , (4.6) 

where ip(tp, v) is the wavefunction evolved in the internal- 
time parameter v, and expressed in the representation 
where acts by multiplication: §i\)(tp,v) = tpip{tp,v). 
More precisely, if the expression on the right of (|4.2j) can 
be interpreted as a conditional probability, then it must 
be possible to interpret (|4.6[) in the same way. Up to 
this point, we have maintained a notational distinction 
between quantities that prima facie have quite different 
meanings, namely {v, tp}, which are internal-time param- 
eters, and {v, tp}, which are configuration-space coordi- 
nates. However, if the internal-time formalism is to be 
interpreted as expressing an observable correlation be- 
tween, in this case, the volume and the scalar field, then 
it must be possible to use these variables interchangeably, 
and we will now do so, until further notice. In particu- 
lar, we would like to identify the probability density (|4.6[) 
with the one displayed in (|4.5p . 

In the light of the arguments summarized at the be- 
ginning of this section, one might expect the conditional 
probability interpretation to present some difficulty. Let 
us, indeed, attempt to construct the required joint prob- 
ability V(y,tp) from the conditional probabilities V(y\tp) 
and V(tp\v), which can be calculated once a wavefunction 
is specified. Denote by f(tp) the denominator in (|4.4[) and 
by g(y) the denominator in (|4.5I) . Then we have 



f{tp)V(y\tp)=V(y,tp) = g{v)2{y\v). 



(4.7) 



Up to constants fo := /(0) and g := ff(0), which are 
fixed by normalization, the unknown functions are deter- 
mined by 

m = 9 4^ (4.8) 



V{Q\tp) 
/o7>HO) 



(4.9) 



WW) 

We see that the joint probability density is well defined 
(the first and last expressions in (|4.7[) are consistent) only 
if 



R 



i. 



(4.10) 



g*V{v\tp)V{Q\v)V{tp\Q) 

In the following, we will calculate the ratio R for a specific 
state, and find that it is not equal to 1. This counterex- 
ample is sufficient to demonstrate that R is not equal to 
1 in general, but we will also argue that our example is 
not especially atypical. 
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B. Inconsistent probabilities 

We begin by constructing the wavefunction ip{u, <f) 
that appears in (|4.1[) and (|4. 2[) . To simplify matters, 
we consider a state for which tp_ — 0, a condition that is 
preserved by evolution with respect to both H v and H u . 
That is, we focus on the sector that corresponds classi- 
cally to an expanding universe, and we will drop the sub- 
script + that indicates this explicitly. The representation 
in which the volume operator V acts by multiplication is 
obtained by the Hankel transform 



de D JotAe^V/V^VM, (4.11) 



where Jo is the Bessel function of the first kind. Bessel's 
equation (zd z zd z + z 2 )J (z) = implies that the opera- 
tor V defined in (|3.16p does indeed act on this function 
by multiplication, if ^(cd, <p) = e leD(p / h ip(ei)) belongs to 
the domain on which V is symmetric. Provided that 
ip(er>,ip) possesses an invertible Hankel transform, it is 
straightforward to verify that this transform preserves 
the normalization, J °° f)\ 2 dv = J °° |V'(eD)| 2 deD. 
We consider normalized states of the form 



V(e D ) = [(2/3) 2 " +1 /(2n)! 



,1/2 



-/3cd 



(4.12) 



where n is a positive integer and /3 a real, positive con- 
stant, which satisfy both of these requirements. Our 
main motivation for this choice is that the Hankel trans- 
form can be calculated analytically. In these states, the 
mean energy of the dust is £b = (n+^)/3 _1 and its disper- 
sion is AeD = (2n + l) _1//2 eD, so for large values of n the 
energy distribution becomes quite sharply peaked. How- 
ever, the analytic expressions for tp{v, ip) become very 
cumbersome for large n, so we focus on the case n = 1. 
In that case, we find 



V{v\ip) = (A 2 //3 7 3 )( 7 - v + \v 2 )^'\ (4.13) 



with 



v := A 2 ^/2/3, 7 := 1 + [ip/phf 



It is easy to check that J °° V{v\tp)&v = 1. 

Evolution in the internal time v 1 ! 2 is generated by £f„, 
whose action on the a = +1 subspace, and in the e D rep- 
resentation is multiplication by (^cd/SttG) 1 ^ 2 — AfieJ/ 2 . 
The representation in which <fr acts by multiplication is 
obtained by Fourier transformation: 

/>OC 

4>(ip, v) = (2Trh)-^ 2 / de D e i£Dv/ri e- iAe D ""Vfe). 
Jo 

(4.14) 

Again, one may easily check that <$>ip(tp,v) = (pip(tp,i/), 

provided that VK £ d) = when cd is or oo, which is the 
condition for the operator (j3.19[) to be symmetric. We 
note in passing that classically, when the constraint is 



satisfied, €d = implies p v — 0. On this hypersurface, 
the volume is constant and the Dirac observable (I2.15[) 
is not well defined. The need to restrict attention here 
to wavefunctions that vanish at 6d = is therefore not 
surprising. 

Using the wavefunction (|4.12p with n = 1, we obtain 



2/? 3 |X(P)| 2 

irh {p 2 + ip 2 /h 2 ) 2: 



(4.15) 



where p := \ 2 is/(f3 — iip/h) and 



1 



X (p) := -^8-2p + iV^(p-6)e-"/ 4 



1-erf 



2 

(4.16) 

We do not know (and neither does any computer-algebra 
package available to us) how to compute the integral 
V((p\v)d<p analytically, but numerical evaluation 
yields the value 1 for randomly selected values of v. 

With these probabilities in hand, we find the ratio 
(|4~TU)) to be 



R = 



7 e- p | X (p)| 2 (l-^ + ^ 2 



-*/7| X (2P)|2( 7 -P+I^) 



(4.17) 



While R is equal to 1 by construction when ip = or 
v = 0, it is not equal to 1 elsewhere. Consequently, the 
functions V(y\<p) and V(ip\v) cannot consistently be in- 
terpreted as conditional probability densities arising from 
some underlying joint probability distribution. Certainly, 
this conclusion is based on a single counterexample, but it 
seems that our sample wavefunction has no special patho- 
logical feature, and the inconsistency is very likely to be 
generic. This will become a little clearer on examination 
of the kinds of evolution that are associated with the two 
internal times ip and v. 



C. Singularity resolution 

It is of considerable interest to see exactly what is 
implied by the two probability distributions. To this 
end, we evaluate them using a wavefunction of the form 
(I4.12j) with n = 4. The resulting analytic expressions 
are lengthy and unilluminating, but the somewhat more 
sharply peaked energy distribution leads to probability 
densities whose nature is more readily apparent to the 
eye. 

Figure [1] shows the probability density V(y\ip) for the 
volume, evolved with the internal-time parameter if. (We 
revert to the notation of (|4.2I) . since the conditional- 
probability notation has proved to be inappropriate.) 
Clearly, the classical singularity at v = is resolved, 
in the sense that it has been replaced by a bounce in the 
quantum theory. (Since the generator H v of evolution 
in ip is independent of a, the probability density for the 
corresponding state in the classically contracting sector 
a = — 1 is exactly the same.) This agrees qualitatively 
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FIG. 1. The probability density for the volume v evolved in 
the internal-time parameter ip. The wavefunction is of the 
form (14. 1 2 p with n = 4, and inessential constants have the 
values h = A = P = 1. 



with the results of Amemiya and Koike 23] , who studied a 
similar model, adopting the Brown-Kuchaf scalar held as 
an internal time, though the quantization schemes they 
considered differ in detail from ours. In loop-quantum- 
gravity- insp ired treatments, such as those described in 
[2l \M 124, the singularity is also seen to be resolved, 
but the mechanism app ears to be different. In partic- 
ular, it is found in [21| (where the matter content is a 
conventional massless scalar field, which is also used as 
internal time) that the minimum volume at the bounce 
corresponds to a density p crit = 3 / (16ir 2 ( 3 G 2 h) , C being 
the Barbero-Immirzi parameter, independent of the de- 
tails of the quantum state. In the present treatment, by 
contrast, the minimum value of (V(tp)} is just (V), which 
is given by 2/3/A 2 in any state of the form (|4.12p . (That 
the bounce occurs at tp — is a consequence of choos- 
ing f3 to be real in (|4.12[) ; it apparent from (j4. 1 1[) that 
an imaginary part of j3 simply shifts ip by a constant.) 



This gives a maximum density p^ 



; ,/(2n+l). 



Physically, a universal critical density of the order of the 
Planck density seems more reasonable, so it is worth em- 
phasizing that our purpose here is to explore the nature 
of time evolution, not to construct an optimal model of 
cosmology, or an optimal quantization scheme. 

Figure [2] shows the probability density V(<p; v) for the 
scalar held, evolved with the internal-time parameter v. 
It clearly indicates a universe expanding from the initial 
singularity: the scalar field, which classically increases 
linearly with proper time, here increases with the volume. 
In the same sense, a state in the a = — 1 sector would be 
seen to contract towards a final singularity. (Specifically, 
replacing ej-/ 2 with — e^ 2 in (|4.14[) leads to the mirror 
image of figure [2j) From this perspective, the two sec- 
tors of the quantum theory reproduce the expanding and 
contracting regions of the classical phase space. A similar 
result is found in 2l| for a Wheeler-de Witt quantization 




FIG. 2. The probability density for the scalar field <p evolved 
in the internal-time parameter v. The quantum state is ex- 
actly the same as in FIG. [1] 



of the model whose matter content is a massless scalar 
field, and more recently in [2o] j. where the same model is 
treated in a consistent-histories approach. 

It is obvious from these figures that the two proba- 
bility densities they depict cannot arise as conditional 
probabilities from an underlying joint probability den- 
sity. We emphasize that they are computed from exactly 
the same quantum state; the difference arises solely from 
the two different notions of internal time used to con- 
struct the families of Dirac observables V((p) and <&{v) 
and their associated Schrodinger-picture wavefunctions. 
It seems clear that these two notions of of internal time 
'see' the expanding and contracting sectors of the theory 
very differently, and that the difference is unlikely to be 
attributable to specific features of our chosen states. 



D. Difficulties with internal time 

The probability densities shown in figures [T] and [2] un- 
derline in a rather striking way (which we did not antic- 
ipate at the outset of this investigation) the difficulties 
of interpretation of the internal-time formalism discussed 
in section IIV Al It seems to us that the two notions of 
internal time must stand or fall together, since they are 
merely two different implementations of the same algo- 
rithm. That is to say, since the two families of Dirac 
observables are constructed in the same way, one cannot 
consistently maintain that V(ip) represents the volume 
when the scalar field has the value ip without accept- 
ing that, by the same token, $(f) must represent the 
scalar field when the volume is v. Since the two associ- 
ated probability distributions, which arise in exactly the 
same quantum state are seen to be in gross conflict, we 



9 



conclude that both viewpoints cannot simultaneously be 
correct, and therefore that neither of them is correct. As 
argued above, the underlying reason for this is that nei- 
ther if nor v can be interpreted as a value assumed by 
any physically observable quantity. 

Some might wish to claim that the conflict should be 
resolved in favour of if as a preferred internal time, be- 
cause the model can be 'deparametrized' in this vari- 
able: that is, the constraint (|2.9[) can be solved to obtain 
V>t> = ~Q( v iPv), where the function is independent of 
</>, and serves as the generator of evolution in the internal 
time if. We do not think that this is a good argument in 
itself, but defer discussion of this point to section [VlJ 

A pressing question raised by the above results is, of 
course, does the quantum theory resolve the classical sin- 
gularity or not? More generally, since in this example 
the answer seems to depend on an arbitrary choice of the 
variable that is to serve as internal time, do we have a 
reliable way of deciding whether or not the singularity 
is resolved in the context of any cosmological model and 
quantization scheme? A possible answer is that quan- 
tum mechanics is actually ambivalent on this question. 
It could be, for example, that given some definite quan- 
tum state, different classes of observer will unavoidably 
disagree on whether that state involves a singularity or 
not. We proposed in [18j . following earlier work in [26| . a 
variant of the idea of relational time that seems to avoid 
the difficulties of interpretation we have stressed up to 
now, by introducing a 'test clock', which provides a pre- 
ferred notion of time evolution from the point of view of 
a specific observer. As described in the following section, 
this preferred description, from the point of view of a 
comoving observer internal to the model universe, effec- 
tively coincides with that furnished by ip as an internal 
time. 



V. DUST-FILLED COSMOLOGY WITH A TEST 
CLOCK 

A. Quantization of an extended model 

As explained at greater length in [l8j], we supplement 
the model studied in previous sections with a rough-and- 
ready description of a small clock, which we consider to 
be internal to a specific observer, and thus localized on 
that observer's worldline. In a complete description, the 
coordinate functions, say x fJ '(s), that locate this world- 
line should appear as extra phase-space coordinates (as 
should a complete set of metric and matter fields, g^(x) 
and (f> a (x)), but in the spirit of simplified cosmological 
models of the sort considered here, we assume that these 
degrees of freedom can be neglected. If the observer is 
comoving, then the proper time along the worldline is 
the same t as appears in (|2.5I) . From the set of vari- 
ables qi that describe the internal workings of the clock, 
we suppose that a function r({q}) can be constructed, 
which constitutes the reading of the clock. Crucially, 



however, we take the view that this reading, being in- 
ternal to the observer in question, is in principle inac- 
cessible to that observer, and for that reason need not 
feature in the physical phase space that describes the 
universe from that observer's point of view. Rather, it 
provides the context for observations of other physical 
quantities to be made. Solution of the equations of mo- 
tion for the variables q~i({q},t) yields the clock reading 

r ({ a }it) '■— r ({<z({<?}^)}) at tne proper time t. This 
reading need not be linear in t, but, if the clock is fit for 
our purpose, there must be, classically, a unique func- 
tion to({q}), such that r({q} 7 t ) — 0. That is, t ({q}) 
is the proper time at which the clock reads 0. It is not 
hard to show that {to, h} = —I, where Nh is the clock's 
Hamiltonian. Because the system described by h is (ide- 
ally) localized on a single worldline, h and to are inde- 
pendent of the cosmological variables (v,p v , <\>,p^). The 
total Hamiltonian constraint is 

C := Co + h = C grav + p 4> + h (5.1) 

and, as shown in [l8j], the quantity 

V(t) := v(v,p v ,(j>,p <t ,\t +t) (5.2) 

is, for each r, a Dirac observable, {V(t),C} — 0. So, of 
course, are $(t), etc., defined in the same way. Thus, 
we obtain a new set of evolving constants of the mo- 
tion, which can be interpreted classically as the volume, 
etc., when a proper time r has elapsed along the ob- 
server's worldline since the clock read 0. The second 
crucial feature of this construction is that r is not to be 
interpreted as a reading obtained from observation of a 
physical clock. From the point of view of the dynamical 
system governed by the constraint (15.11) , it is an external, 
'heraclitian' time, in the sense of Unruh and Wald 27]. 

We quantize this enlarged model using essentially the 
same scheme as in section lllTl adopting a kinematical rep- 
resentation in which i, p^ and h act by multiplication on 
wavefunctions ^(v, er>, e c ), where e c is the energy of the 
clock. The kinematical operator representing the fidu- 
cial proper time to is then to = —ihd/de c . Since we are 
no longer using the notion of evolution in the internal- 
time parameter z/, no operator of interest distinguishes 
the two sectors a = ±1, and we will focus on a single 
sector, writing a solution to the constraint equation as 

e D , e c ) = ^(e, 6 D )C (Ae 1 / 2 U 1 / 2 ), (5.3) 

where e := £d + £c is the total energy, and Co is any 
Bessel function of order 0. The physical Hilbert space is 
now Hphys = £ 2 (K+, deDde c ), with the inner product 

(V'l^phys 

= de D / de c f/ii(e D + e c ,e D )-02(eD + £c,£d)- 
Jo Jo 

= / de / de D ^i(e, e D )tp 2 (e, e D ). (5.4) 
Jo Jo 
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The intention is that e c should be small compared with 
£d in the same sense that the biological system respon- 
sible for an astronomer's ciracdian rhythm, say, is small 
compared with the energy content of the visible universe, 
and we implement this by restricting attention to states 
such that \ip\ 2 is very small unless e c < cd- The result- 
ing limitations on the resolution with which r-dependent 
observables might be determined are examined in some 
detail in [l8[ for a universe whose matter content is a 
massless scalar field, and we will not repeat the analysis 
for the present model. 

Acting in H p hy S , we find families of Dirac observables 



V(t) = V- 8{h\)- 2 YT - A(hX)- 2 C gia 

<S(Y) = l> + T, 

where, acting on ip(e, £d), 

V = -(4/A 2 )<9 e e<9 £ 



(5.5) 
(5.6) 



Cgrav — 



$ = ihd e 



(5.7) 
(5.8) 
(5.9) 
(5.10) 

Note, in particular, that in (I5.10[) . the derivative is with 
respect to £d keeping the total energy e fixed. Treating eD 
and e c as independent variables, we have $ = <\> + to = 
ih(d eD — <9 £c ). Consequently, the observables V(t) and 
$(t) now commute for every r and, according to the 
usual rules of quantum mechanics, can simultaneously be 
assigned measured values y_ and tp. A domain on which 
all of these operators are symmetric is provided by the 
boundary conditions 



^(e,e) = V(e,0) = S^(e, e D )| e =e D = 0. 



(5.11) 



Of course, e = cd means that the clock energy e c van- 
ishes, and we always take ip(e, cd) to vanish fast enough 
at large arguments to ensure square integrability. 



B. Joint probability density 

Evolution of the Heisenberg-picture operators ([57 
and (|5.6|) is generated by the self-adjoint Hamiltonian 

H T = Cg Iav 

+ P<t> = £d — e, and states satisfying the 
boundary conditions (|5.11l) can be expressed in a repre- 
sentation in which both V and $ act by multiplication, 
through combined Fourier and Hankel transformation. 
Thus, we can now define a bona fide joint probability 
density, which evolves in the usual way with proper time 
t, namely 



ViMsWr) = \ip(v,<f, T )\ 2 , 



(5.12) 



where 



x 



/ de / d£ D e i^v/h e i(e-e D )r/h 

Jo Jo 

(5.13) 



JoCAeVV/'MceD) 



In particular, this probability density has the r- 
independent normalization J Q dvJ_ oo dipV(v,ip;T) = 1. 

To make contact with the probabilities discussed in 
section HVl consider a state that factorizes as 



(5.14) 



where ip c is a wavefunction for the clock and, as above, 
e c = e — £d is the clock's energy. (Recall, though, that 
the clock is not represented by any operator independent 
of V and $ acting in 'Hphys, and r is not a value obtained 
from observation of the clock). The wavefunction (|5.13p 
becomes 

$<JL®r) = ± / dee^MXe 1 ^ 2 ^ 2 )^) 
z Jo 

de c c ie ^ T -^/ r Vc(ec). (5.15) 

The idea of a test clock described previously implies that 
tp(e) should be peaked around a value e of the order of the 
energy content of the visible universe, while ipc( e c) should 
be peaked around a value e c somewhat smaller than the 
mass of an astronomer. Under these circumstances, it is 
an excellent approximation to extend the upper limit of 
the e c integral to infinity, in which case the wavefunction 
(|5.15j) and the probability density (|5.12[) also factorize. 
In fact, we have 



V(v,<p;t)~V(js\<p)V((p;t), 



(5.16) 



where V(v\tp) coincides with the probability density de- 
fined in (|4.2|) . except that we do not now need to distin- 
guish the sectors a = ±1. Here, however, the arguments 
y_ and <p genuinely stand for values obtained from obser- 
vation of the quantities represented by the commuting 
operators V and $, and V(v\tp) is a genuine conditional 
probability. The function 



1 

2~Trh 



= {T-<p)/fr 



V>c(e c ) 



(5.17) 



is just the probability density for obtaining the value tp 
from a measurement of the scalar field at proper time r. 
Given a suitable choice of ip c {e c ), it will be sharply peaked 
on a trajectory of the form ip = $o + t, consistent with 
the solution (|5.6p of the Heisenberg equation of motion. 



C. Timeless interpretation 

We have emphasized that the parameter r, which la- 
bels the families of Dirac observables V(r) and $(t) is 
an unobservable, external time parameter. Classically, it 
coincides with the arc length of an observer's worldlinc, 
on which the fiducial event, that the unobserved clock 
reads 0, serves to define an origin. Naturally, the time- 
dependent observations recorded by an astronomer who 
wishes to test a theoretical expression such as (|5.12l) will 
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not refer either to r or to the internal, unobserved clock 
that we have pictured, for the sake of argument, as a 
biological clock. Instead, they will refer to the observed 
readings of some time-keeping device, which we will call 
(taking a cautious view of the generosity of the relevant 
funding agency) a 'wristwatch'. 

From an operational point of view (more or less in the 
sense of Bridgman[28|) substantive physics is contained 
only in correlations between measured values of observ- 
able quantities. The idea that a solution to the problem 
of time in constrained systems is to be sought in such 
correlations is explicitly developed in, for example, [29i — 
l3lj . and is implicit in much of the recent literature. We 
have argued that the required correlations are not di- 
rectly provided by a wavefunction such as (|4.ip , but in 
the model considered here, it is straightforward to see 
that this wavefunction does indirectly lead to to an es- 
timate of the desired correlation, to an approximation 
that might be extremely good. In principle, we would 
like to study correlations between observed values of the 
volume, the scalar field and the astronomer's wristwatch. 
Since the wristwatch is another small clock, which clas- 
sically follows, for practical purposes, the same worldlinc 
as the unobserved test clock, we could incorporate it into 
our model by adding its energy h w to the total constraint: 

Ctotal := Cgray + + fl w + h. (5.18) 

This presents no technical difficulty, but since, for this 
particular model, p^ and h w appear in the same way 
in Ctotal, it is clear that we can usefully economize on 
clocks by deleting h w and treating the scalar field as 
the astronomer's timekeeping device. (Here, we profit 
from the enormous simplification that results from the 
assumption of homogeneity. In a more general setting, we 
envisage that Dirac observables depending on t can be 
constructed only from fields in the observer's immediate 
locality. Technical difficulties aside, this is no significant 
limitation, since the local fields include, for example, the 
CMBR photons entering the astronomer's telescope.) 

Suppose, then, that simultaneous measurements of vol- 
ume and scalar field are performed at a sequence of 
proper times Tj, distributed according to some density 
function T)(t), with J_ r7(r)dr = 1. (t](t) may have 
support only in a subinterval of (—00, 00) if, for example, 
the observer's worldline terminates at a singularity.) The 
function 

/oo 
V{v,^r)r,{T)dr (5.19) 
-00 

is a correctly normalized joint probability density, which 
furnishes a timeless description of the correlation between 
these measured quantities In the approximation that 



2 This argument takes no account of any 'collapse of the wave- 
function' occasioned by performance of the measurements. Con- 



the joint probability density factorizes as in (|5.16[) . we 
find 

(5-20) 

where V((f) = J^° QO 'P( ( P]T)r](T)dT is the probability for 
finding the measurement of the scalar field on a randomly 
selected occasion to have yielded the value if. For the 
state studied in section lTVCj V(v\<p) is precisely the func- 
tion depicted in figure[T] but the introduction of an unob- 
served test clock now allows us to interpret this function 
as a bona fide conditional probability. In that indirect 
and approximate sense, we identify <p as a preferred in- 
ternal time, and conclude that the singularity is resolved 
in this quantum theory — or at least that the observer 
whose unobserved clock is represented by h will discover 
herself to be living in a bouncing universe. 

In part, this conclusion is specific to our somewhat 
artificial model, in which the matter content of the uni- 
verse is provided by the Brown-Kuchaf scalar field. More 
generally, consider a Hamiltonian constraint of the form 

C := C grav+mat t + Pw + h = 0. (5.21) 

In this expression, C grav +matt is the contribution of met- 
ric and matter fields, h is the Hamiltonian of an unob- 
served test clock, and p w is the momentum conjugate to 
the pointer reading r w of a wristwatch following essen- 
tially the same worldline as the test clock. By taking the 
Hamiltonian of the wristwatch to be just p w , we model a 
time-keeping device that is manufactured so as to supply 
a linear measure of the proper time along its worldlinc. 
Following the same steps as before, we expect to obtain 
a physical Hilbert space of functions ip(~Y,e,e w ), where 
e = e w + e c is the total energy of the wristwatch and 
the unobserved clock, while 7 collectively denotes the re- 
maining metric and matter variables. Operators on 'Hphys 
representing Dirac observables include 

Cgrav+matt = — ( (5.22) 

r w = ih^-, (5.23) 
de w 

as in (|5.9p and (|5.10p . Taking a wave function that fac- 
torizes as -0(7, e, e w ) = ^(7, e)-0 c (e c ), we will obtain a 



sequently, in the spirit of the Copenhagen interpretation, the 
sequence of measurements should strictly be understood as be- 
ing performed on an ensemble of identically prepared universes. 
We find this deeply unsatisfactory, but the problem concerns the 
interpretation of quantum mechanics, especially as applied to the 
universe as a whole, about which we have nothing useful to say, 
rather than the interpretation of time, about which we believe 
we do have something useful to say. We adopt a Copenhagen-like 
point of view, not because we find it convincing, but because we 
find rival interpretations to be more cumbersome without being 
any more convincing. While recognizing that many physicists 
would disagree with this, we do not think that the present paper 
would be enhanced by a digression on this issue. 
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Schrodinger-picture wavefunction analog ous to (HHSJl, of 
the form 

/>00 

^(7,r w ;r)= / dee^w/^^,) 
Jo 

x f de c e k ^-rw)/S/, c ( ec ). ( 5 . 24 ) 
Jo 

From this, we obtain a joint probability P(7, £ w ;t) = 

l^(T)£wi r )l 2 ano -i ^ we w ishj the timeless version 
■p(7,r w ) = J 7^(7, r w ; T)7j(r)dT, which describe correla- 
tions between observed values of the cosmological quanti- 
ties 7 and the wristwatch pointer r w . However, because 
the energy e = e w + e c is now the total energy of two 
small clocks, the wide separation of energy scales which 
led to the factorization in (I5.16[) is no longer present, and 
the cosmological field ip is replaced by the wristwatch 
pointer r w . Generically, therefore, we do not automat- 
ically recover a conditional-probability interpretation of 
any internal time chosen from among the cosmological 
variables 7. 



VI. DISCUSSION 

We argued in [l8| that the notion of internal time, es- 
pecially as commonly implemented in models of quantum 
cosmology, is unsatisfactory, for two reasons. First, it of- 
fers no account of the passage of time as it is ordinarily 
conceived. Ordinarily, there seems to be a clear sense 
in which a well-constructed clock reads '10s' seven sec- 
onds after it read '3s', and this does not appear merely 
to result from a conspiracy amongst the manufacturers 
of time-pieces. That is, a time-keeping device can be said 
to work accurately (or not), because there is a time to 
be kept, not merely because its readings tend to agree 
(or not) with those of other devices of the same sort. 
In general relativity, this notion of time is provided, for 
some specific observer, by the proper time that elapses 
along that observer's worldline, but it is not recovered in 
a treatment that describes evolution by correlating, for 
example, the volume of a spatial region with the value of 
a scalar field. Second, as summarized in section lTV A[ the 
meaning of a parameter that serves as an 'internal time' 
is unclear; in particular, it cannot be construed as a value 
obtained from the observation of any physical quantity 
that is to be regarded as a clock, because no such quan- 
tity is represented by any operator acting in the physical 
Hilbert space H p hys- 

The textbook interpretations of quantum mechanics 
are not, of course, designed to deal with dynamics gen- 
erated by a Hamiltonian constraint, and some modifica- 
tion might well be needed to deal with that situation. 
In particular, one should perhaps reconsider the usual 
statements that observable quantities are represented by 
symmetric operators on Hphys: an d that the values of two 
such observables can be determined simultaneously only 
if their corresponding operators commute. Concretely, 



for models of the kind considered here, it is certainly 
possible to compute a wavefunction 1(1(1/, (p), which os- 
tensibly evolves with respect to internal time ip, and one 
would like to be able to interpret it as expressing a cor- 
relation between observed values of the volume and the 
scalar field. In view of the ^-independent normalization 
J \ip(v, ip)\ 2 du = 1, the quantity \ip{v, ip)\ 2 cannot be in- 
terpreted as a joint probability density, but it might be 
possible to regard it as a conditional probability density 
for the volume, given that the scalar field has the value 
ip. If so, then (i) there must be some underlying joint 
probability density, not directly supplied by the usual in- 
terpretation of the wavefunction, and (ii) there must exist 
a complementary conditional probability density for the 
scalar field, given some value for the volume. 

To investigate this possibility, we have studied a sim- 
ple model and quantization scheme, whose main virtue is 
that it admits the explicit construction within the same 
quantization scheme of two complementary internal-time 
evolutions. We found, though, that the two comple- 
mentary probability distributions, computed for the same 
quantum state are inconsistent. In fact, these two notions 
of internal time lead to qualitatively different physical in- 
terpretations of the same quantum state: with the scalar 
field as internal time, the universe appears to bounce 
(figure [I} , whereas with the volume as internal time it 
appears to expand from an initial singularity (figure [5]) . 
In view of this inconsistency, we conclude that neither 
of these two notions of internal time can be correct. If 
the operator V(ip) really represents the volume when the 
scalar field is tp, then it must also be true that &(v) is 
the scalar field when the volume is and that cannot be 
so. 

We now wish to conclude further that the notion of 
internal time, as implemented by the evolving constant 
of the motion construction, is not tenable in general. To 
be sure, we have discovered a serious inconsistency only 
in the context of one particular model and one particular 
quantization scheme for which the required computations 
are straightforward, but if the idea were sound in general, 
then it ought to apply to this example. 

As described in section [V] a variant of the relational- 
time construction proposed in [l|| [2(| is capable of cir- 
cumventing these difficulties. We augmented the model 
by including an idealized description of a small test clock, 
which we imagine to be internal to a specific observer, in 
this case a co-moving observer. The time parameter r 
that labels families of Dirac observables is, classically, at 
least, the proper time that elapses along the observer's 
worldline. Its value is not to be regarded as obtained from 
observation of any physical clock and, from the point of 
view of the observer in question, plays the same role as 
the external time in textbook Newtonian or quantum 
mechanics. The test clock itself we take to be unobserv- 
able in principle by the observer in question, by virtue 
of being internal to that observer. Using this construc- 
tion, we could obtain bona fide joint probabilities, which 
describe genuine correlations between cosmological ob- 
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servables and whatever time-keeping the observer might 
use in the course of recording observations. 

In the case that the observer uses the cosmological 
scalar held as a (large) clock, the resulting conditional 
probability density for the volume essentially coincides 
with the evolving probability density obtained using the 
scalar field as an internal time. This is a special fea- 
ture of the Brown-Kuchaf field, whose contribution to the 
Hamiltonian constraint (I2.9[) is just its canonical momen- 
tum. For the same reason, this model is deparametrizable 
in the variable 4>: in effect, by substituting p^ — > — i/icL, 
one converts the constraint equation into a Schrodinger- 
like equation that governs evolution in the internal time 
(p. We do not think, however, that deparametrizabil- 
ity is in itself sufficient to identify a preferred internal 
time, though it does, of course, make the implementation 
straightforward. In general, a model is deparametriz- 
able in a variable <fi if the constraint C(<p,p^, w) = can 
be solved to obtain p^ — —Q(uj), where is indepen- 
dent of </>, and u) denotes the remaining canonical vari- 
ables. Again, one can substitute p^, —¥ —ihd v to obtain 
a Schrodinger-like equation, but, again, the parameter ip 
cannot be interpreted as the observed value of some phys- 
ical quantity, and it is hard to see what other meaning 
it might have. Moreover, this strategy suffers from what 
Kuchaf gj calls the 'Hilbert space problem'. That is, the 
Schrodinger-like evolution can be implemented only if 
the inner product on "H p hy S is chosen so as to make 
self-adjoint. The same inner product will not necessar- 
ily confer self-adjointness on the generators of evolution 
with respect to other candidates for an internal time, and 
it seems somewhat unreasonable that the inner product, 
and hence the quantum theory as a whole, should de- 
pend on this arbitrary choice of an internal time. In the 
present example, <j) has a preferred status not because of 
deparametrizability as such, but rather because its con- 
tribution to the constraint is linear in with a constant 
coefficient. In fact, it is worth noting that the canoni- 
cal transformation u = v 1 / 2 and p u = 2v 1 ^ 2 p v brings the 
constraint (|2.9I) to the form Co = — (Xh)~ 2 p 2 l +p l j > = 0, so 



that this model is also deparametrizable in the variable u. 
(Some authors contrast a 'relativistic' deparametrization 
in u, meaning that Co is quadratic in p u , with a 'non- 
relativistic' deparametrization in <j), meaning that Co is 
linear in p$.) This simple form of Co, depending only 
on momenta, explains why it is straightforward, in this 
model, to obtain the two notions of evolution with re- 
spect to internal time (|2.23[) - (|2.26[) . with generators that 
can simultaneously be made self-adjoint in the quantum 
theory. 

We believe that the idea of a test clock, internal to 
some specific observer, as implemented here provides an 
improved notion of time evolution, but it would be at 
best premature to suggest that it yields a definitive so- 
lution to the problem of time in general. Among the 
limitations of the proposal as we have so far described 
it are the following. (1) While the construction appears 
to be successful in spatially homogeneous models with a 
single constraint, it does not follow that a similar con- 
struction will work in more general spacetimes. (2) In 
particular, we have bypassed any explicit description of 
the observer's worldline by considering only a comoving 
observer in a Friedmann-Robertson- Walker universe. (3) 
The time parameter r survives the passage to quantum 
mechanics as a c-number parameter, but it is not obvious 
that this parameter can be unambiguously described in 
the quantum theory as the arc length of a worldline. (4) 
In the context of simple cosmological models, at least, it 
seems to be an inevitable consequence of the constraint 
that some object to which one is inclined to attribute 
a real physical existence turns out to be unobservable. 
We think it is plausible that a clock which is internal to 
some observer should turn out to be unobservable in a 
description of the universe 'from that observer's point of 
view'. However, this plausible form of words is not di- 
rectly mandated by the formalism, and some other way 
of understanding the unobservability of whatever quan- 
tity is eliminated by solution of the constraint may turn 
out to be better founded. We plan to address these issues 
in future work. 
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